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ABSTRACT 

The possibihty that long tidal tails formed during compact object mergers may power optical tran- 
sients through the decay of freshly synthesized r-process material is investigated. Precise modeling 
of the merger dynamics allows for a realistic determination of the thermodynamic conditions in the 
ejected debris. The results of hydrodynamic and full nuclear network calculations are combined to 
calculate the resultant r-process abundances and the heating of the material by their decays. The sub- 
sequent homologous structure is mapped into a radiative transfer code to synthesize emergent model 
light curves and determine how their properties (variability and color evolution) depend on the mass 
ratio and orientation of the merging binary. The radiation emanating from the ejected debris, though 
less spectacular than a typical supernova, should be observable in transient surveys and we estimate 
the associated detection rates. The case for (or against) compact object mergers as the progenitors 
of short gamma-ray bursts can be tested if such electromagnetic transients are detected (or not) in 
coincidence with some bursts, although they may be obscured by on-axis afterglows. 
Subject headings: nuclear reactions, nucleosynthesis, abundances — black hole physics — radiative 
transfer — stars: neutron — hydrodynamics — gamma-ray burst: general 



L INTRODUCTION 

Merging compact binaries are the primary candidate 
for dire ct detection of grav itational waves (GWs) by 
LIGO dAbbott et"ar]|2QQ8D , and are thought to be 
the progenitors ot short gam ma-ray bursts (SGRBs)(c.f. 
[Lee fc Ra mirez-RuTz 2007; Rosswog 2007; Nakar 2007j 
|Gehrels"t al.^^2009) . Current observational limits indi- 
cate that any supernova-like event accompanying SGRBs 
would have to be over 50 times fainter than normal Type 
la SNe or Type Ic supernovae (SNe), ( [Hjorth et al |20Q5 



[Bloom et al 2006). These limits strongly constrain pro- 
genitor models fo r SGRBs ( Lee fc Rami rez-Ruiz 2007; 
Perley et al1|2009[ [Kocevski et aT.||2UlU|) . Unless SGRBs 
are eventually found to be accompanied by telltale opti- 
cal signatures like th e supernovae of long-duration GRBs 
(Hjorth et al. 2003), the only definitive understanding 
of the progenitors will come from possible associations 
to direct gra vitational or neutrino signals ( [Rosswog fc| 
Liebendorferl[2003) . 



I'he merger of compact objects does not necessarily 
imply the absence of optical or other long-wavelength 
phenomena after coalescence. Neutron-rich material may 
be dynamically ejected during a NS-NS (Neutron Star) 
or a NS-BH (Black Hole) merger. Material dynami- 
cally stripped from a star is violently ejected by tidal 
torques through the outer Lagrange point, removing en- 
ergy and angular momentum and forming a large tidal 
tail. Its subsequent decompression may synthesiz e ra- 
dioactive elements through the r-process (Lattimer fc 
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'Schramm' '1976J [Freiburghaus et al. 1999), whos e ra- 
dioactive decay could power an optica l transient (Li fc 
iPaczyh ski 1998] [Metzger et al.|[2QlQbl ). For doubled 
binaries there are one or two such structures, depending 
on the m ass r atio an d equation of state (EoS) of the NSs 



( Oechslin et al^007 ) . Clearly, only one tail (or no tail) is 
formed in a BH-NS merger, and these are typically a few 
thousand kilometers in size by the end of the disruption 
event and some of the fiuid (as much as a few hundredths 
of a solar mass) is often gravitationally unbound. 

The most efficient conversion of radioactive energy 
to radiation is provided by those isotopes with a de- 
cay timescale comparable to the radiative diffusion time 
through the ejecta. In reality, there are likely to be a 
large numb er of nuclides with a very broad range of decay 
timescales ( |Li fc Paczyhski||1998[ [Metzger et al.[|2010b[ ). 
Current observational limits thus place interesting con- 
straints on the nuclear evolution of this material, as well 
as the total mass ejected in these events. 

In this Letter^ we present the results of multi- 
dimensional hydrodynamics simulations of neutron star 
mergers for various binary mass ratios. The ejected ma- 
terial is then post-processed with a full nuclear network 
and a radiation transport code. Our goal here is to in- 
vestigate how radioactive r-process elements synthesized 
in the tidal tails could power electromagnetic transients 
and how their properties may depend on the mass ra- 
tio of the merging binary. We also discuss prospects for 
detection and implications for SGRBs. 

2. TIDAL TAIL EVOLUTION 
2.1. Ejection 

To calculate the ejected mass and structure of the tidal 
tails, we employ a three dimensional smoothed parti- 
cle hydrodynamics (SPH) method to compute the evo- 
lution of t h e compact objec t merger (Lee fc Ramirez- 
Ruiz 2007 Lee et al. [2010 ). Due to its Lagrangian 
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Figure 1. Homologous structure of ejected material in the four models considered in this work (see table [T]) viewed face on (top) and 
edge on (bottom). The density contours are not the same between models, we have instead chosen to plot Tne models at the same size 
scale. The peak density contour (blue) is twenty times the minimum density contour (red). The evolution from one to two tidal tails in 
the NS-NS mergers as the mass ratio approaches unity is clearly visible. 



nature, SPH is perfectly suited to follow tidal disrup- 
tion processes during which the corresponding geome- 
try, densities, and time scales are chanfflng violently (c.f . 
Rasio fc Shapiro||1994[ [Lee||2QQQ[ [Rosswog et"aL]|2QQ3D . 
I'he equation ot state is a hybrid, where we have com- 
bined the cold Friedman-Pandharipande-Skyrme (FPS) 
nuclear Eo S with an ide a l ther mal compoment as de- 
scribed by Shibata et al. (2005). Once the fluid in the 
tidal tails has become unbound, the full SPH calculations 
are stopped and the particles are followed on ballistic tra- 
jectories in the potential of the central compact remnant 
until their expansion becomes homologous. 

To explore the possible range of outcomes depending 
on the physical parameters of the binary system, we have 
considered four models: three NS-NS mergers with mass 
ratios q = M2/M1 = 1, 0.95, 0.88 and one BH-NS merger 
with q = 0.31. Around O.O5M0 of material is ejected in 
all four models. Relevant parameters resulting from the 
simulations are given in table [l] The density structure of 
the tails is shown in figure [l] at homology. Clearly, there 
is a progression from equal mass tails in the g = 1.0 case 
to almost no secondary tail in the q = 0.88 NS-NS case. 
The relative mass of the two tails is a function of the 
mass ratio and the nature of the merging objects, as well 
as the underlying EoS (which we have not explored in 
this work). It is of interest to determine if different tail 
geometries are distinguishable in the optical signal from 
a compact object merger, as this would provide a signif- 
icant constraint on the nature of the objects involved in 
addition to the expected gravitational wave signal. 

2.2. Nuclear Evolution 

We follow the evolution of the composition of the tails 
using a 6312 isotope nuclear reaction network which ex- 
tends past Uranium. Density histories of particles from 




Figure 2. Top: energy deposition rate from nuclear decay (includ- 
ing neutrino losses) as a function of time for various Lagrangian 
trajectories from the SPH simulations. These are shown for the 
BH-NS merger, but are representative of the NS-NS mergers as 
well. The trajectories are color coded by their density one day 
after the explosion. The gray lines show the heating rate from sin- 
gle reactions that contribute significantly after 0.1 days for a single 
trajectory. Bottom: final abundances as a function of nuclear mass 
for the same trajectories. 
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Table 1 

Characteristics of simulations 



Type 


Mass Ratio 


Primary Mass 


Ejected Mass 


Eject a Velocity 




^peak 


^dec 








(Mq) 


(Mq) 


(c) 


(s) 


(days) 


(days) 


(ergs 1) 


NS-NS 


1.00 


1.4 


0.057 


0.202 


3.2 


0.93 


1.7 


1.51 X lO'^^ 


NS-NS 


0.95 


1.4 


0.047 


0.200 


2.0 


0.93 


1.6 


1.19 X 10^2 


NS-NS 


0.88 


1.5 


0.057 


0.205 


2.1 


1.02 


1.8 


1.44 X 10^2 


BH-NS 


0.31 


5.4 


0.060 


0.248 


3.9 


0.93 


1.7 


1.64 X 10^2 



^ Time to reach 1% deviation from homologous evolution. 
^ Time to reach peak bolometric luminosity after merger. 
^ Light curve decay timescale. 
^ Bolometric Luminosity. 



the SPH simulations are employed. The charged par- 
ticle and neutron capture rat es in the network up to 
At (Z = 85) are taken from [Rauscher fc Thielemann 
(|2QQQ ) . Pas t At, the neutron capture rates ot Fanov 



eraTTpOlO) are employed. The network terminates at 
1U2. Experimental values are taken for nuclear 



masses where available, elsewhere theoretical masses are 
taken from Moller et j il. (2003 ). Neutro n induced fission 
rates are taken fr om | Panov et al. | (|2010,) and the simple 
approximation of l^rankel Metropolis ( 1947[ ) is used to 
calculate s pontaneous fission rates, -b'ission barriers are 
taken from Mamdouh et al.| ([2001). For our fission frag- 
ment distributions, we employ the empirical fits of |WahI| 
([2002). 

I'he nuclear net work is evolved in tinie usin g a variant 
of the code XNet (iHix & Thielemann||1999D. We have 



implemente d the PARDISO sparse matrix solver (Schenk 
et al.[|2008 l), which makes calculations of large networks 



with implicitly coupled fission interactions feasible. The 
energy released by nuclear reactions i s self-con sistently 
added back to the material, similar tO jFreiburghaus et al., 
( [19991 ). 

'I'he material ejected in the tail is from near the surface 
of the neutron star, making it challenging to determine 
the initial electron fraction and entropy of the material 
from our SPH simulations. Therefore, we assume the 
initial electron fraction to be = 0.2 and start the 
calculation at p = 10^ ^g cm~^ and Tg = 1. The ini- 
tial composition is then determined by nuclear statistical 
equilibrium (NSE). As was pointed out by Goriely et al.| 
(^2005), if tidally ejected material is un-shocked, the mi- 
tial temperature of the material is likely much less than 
our assumed temperature. To test the dependence on 
the initial temperature, we have run models going down 
to an initial temperature of Tg = 0.2, and find that the 
nuclear heating rate is not substantially altered. 

The total heating rate and final abundance distribution 
as a function of mass for a n umber of fiuid elements are 
shown in figure [5] Similar to Metzger et al. (2010b), we 
find that the late time heating rate is insensitive to the 
exact initial con ditions and is statisti cal in nature, as 
was predicted by Li & Paczyhski' ('1998') . At one day, the 
top five beta-decays contributing to the heating rate are 
i25Sb, i26sb, 1321^ i27Te, and ^^^Pt. 

2.3. Radiative Transfer 

In order to study the radiative transfer problem, we 
employ the SEDONA 3-D time-dependent LTE Monte 
Carlo radiative transfer code which includes gamma-ray 



transfer and spectropolarization (Kasen et al.'2006). The 
output of the hydrodynamic simulations is mapped into 
the radiative transfer code and the properties of the 
emergent radiation are calculated. A global nuclear heat- 
ing rate based on a fit to the nuclear network calculations 
is employed. We approximately account for neutrino 
losses by assuming 75% of the nuclear network energy 
generation is deposited in the material (Metzge r et al^ 
2010b). Of the energy that is left, we assume 50% is de- 
posited as gamma rays from decays while the other 50% 
is deposited thermally. 

The opacity of r-process material at the relevant densi- 
ties and temperatures is not well known. The main con- 
tribution to the opacity is presumably due to millions of 
atomic lines, which are Doppler broadened by the high 
differential velocities in the eject a. Unfortunately, com- 
plete atomic line lists for these high-Z species are not 
available. Given the uncertainty, we assume in these cal- 
culations a constant gray opacity of = 0.1 cm^ g~^ 
which is characteristic of the line expansion opacit y from 
iron group elements (e.g. |Kasen_ fc Woos leyll2QQ7t. This 
is in contrast to the work of Metzger et al. (2010b), where 
oscillator strengths for pure he were assumed with ion- 
ization potentials from Pb. Considering that neither ap- 
proach will yield accurate spectral information, we feel 
that our simple grey opacity s cheme is as reasonable an 
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2010b). 



approximation as that used in Metzger et al. 

We directly calculated the spatial distribution of 
gamma-ray heating by following the transport of gamma- 
rays and determining the fraction of their energy ther- 
malized by Compton scattering and photoelectric ab- 
sorption. Since the dominant Compton opacity has only 
a weak wavelength dependence, the exact spectrum of 
gamma-ray emission from the radioactive source does not 
strongly affect the results. We therefore simply assumed 
all gamma-rays were emitted at 1 MeV. Before two days, 
the gamma-ray thermalization rate is greater than 80%. 

3. DETAILED PROPERTIES OF THE 
ELECTROMAGNETIC COUNTERPARTS 

In the top panel of figu re [3| the R-band light curves for 
the four models in table Hjare s hown. As would be ex - 
pected from the simple models of Li fc Paczyhskil(1998), 
the peak luminosity correlates with the total ejected mass 
of radioactive elements and the time of the peak scales 
inversely with mass and velocity. Because the total mass 
ejected in these mergers is not very sensitive to the 
nature of the merger cannot easily be determined solely 
from the peak time or luminosity. Additionally, the vari- 
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Figure 3. Top: r-band luminosities for the three models described 
in table^ The solid lines correspond to the q = 1.0 NS-NS merger, 
the dashed lines to the q = 0.95 NS-NS merger, the dashed lines 
to the q = 0.88 NS-NS merger, and the dot-dashed lines to the 
BH-NS merger. The black lines are the r-band luminosity for the 
models (left axis), averaged over all solid angles. The gray lines are 
the g — r color evolution for the models (right axis). The colored 
lines give the luminosity as a function of polar angle (averaged over 
the phi direction) for the q = 1.0 NS-NS merger. Bottom: r-band 
magnitudes for the same models, along with observed SGRB after- 
glow absolute magnitudes (filled squares) and upper limits (down- 
ward pointing triangles) from the compilation of Berger (2010) 
where red shifts were obtained. All of the points correspond to 
separate events. Synthetic SGRB afterglow light curves from van] 
[Eerten &: MacFadyen ( 2011) for a jet energy of 1 x 10^^ ergs are 
also sliown, where tlie solid (dashed) red line is for a jet opening 
angle of 0.2 (0.4) radian. The solid green curve is the optical light 
curve of SN94I, a type Ic supernova. The dashed green curve is for 
SN2001el, a type la supernova. 



ation in peak luminosity with viewing angle within a sin- 
gle model is as large as the variation in the angle averaged 
peak luminosity between models. This is also shown in 
figure 3) This further complicates our ability to distin- 
guish between different mass ratios and progenitor mod- 
els based only on luminosities 

Still, it may be possible to determine if one or two 

^ We also d o not observe the non-s mooth structures seen in the 
light curves Qf |Metzger et al.| ( |2010'b| >, which are due to their use 
of approximate non-grey opacities. 



tails are present based on the color evolution of the 
light curves. In NS-NS mergers which produce two 
tails, the luminosity of the transient will be given by 
the sum of the luminosities of the tails, each of w hich 
can be approximated as a Li & Paczyhski (1998) ex- 
panding sphere. We denote here the heavier and lighter 
tails with the subscript 1 and 2, respectively. |Li & 



|Paczyhski| (1998) find that at late times, the evolution 
of the ettective temperature is given by Teff ^ 4.7 x 

10^ K (M/O.OlM0)lA(c/^)3/4(^^y/^)3/4(^/3 X 10-5)1/4^ 

The ratio of the effective temperatures in the tails at late 
times is then 



^eff,l 



'Ml 

Teffo ' V M2 



1/4 



3/4 



(1) 



Significant variation in color from the single tail case 
is then expected only for a considerable difference be- 
tween velocities. In this case, tail 2 will shift the to- 
tal light blueward if it makes an important contribu- 
tion to the total luminosity at any time. The time 
of peak luminosity for a single tail is given by tm ^ 
Iday (M/O.OlM0)i/2(3v/c)^/^ If the velocity of the sec- 
ond tail is much lower, it contributes more to the total 
luminosity at late times and the light will be bluer com- 
pared to the more massive tail emitting radiation alone. 

In our detailed models, the velocity difference between 
the tails is not significant enough compared to the mass 
difference between the tails in either of the asymmetric 
NS-NS merger models for the tails to be easily discernible 
in their color evolution, as shown in figure [Sj As a re- 
sult, there is no significant distinguishing characteristic 
between eject a geometries in their colors. All four models 
show similar reddening as a function of time. 

Without knowing the detailed line opacity of the R- 
process elements, it is difficult to predict the spectro- 
scopic signatures of NS mergers. We expect that high-Z 
atoms will have a very large number of lines in the opti- 
cal, similar to and perhaps exceeding the number known 
for the iron group elements. Because the ejecta veloci- 
ties in the tidal tails are a factor of ~ 2 — 3, larger than 
ordinary supernovae, the absorption features in the pho- 
tospheric phase (< 15 days) should be very broad. The 
likely outcome is a blending together of the lines into 
a relatively featureless continuum, not unlike the early 
time observati ons of broad line Type Ic supernovae (c.f. 
Galama et aLl[l998 ). 



More detailed information about the ejecta geometry 
could be inferred from spectroscopic observations taken 
at late times (> 15 days) when the ejecta has entered 
the nebular phase and emission features begin to appear. 
Assuming relatively unblended lines (or line complexes) 
are present at these epochs, the spectral features should 
show double peaked emission profiles, at least for equa- 
torial viewing angles. The relative strength of the peaks 
could be used to estimate the relative masses of the two 
tails. 

In brief, the inclusion of realistic ejecta geometries does 
not significantly alter the predictions made by simple 
spherical "tail" models which include realistic nuclear 
physics. This can be attributed to the fact that the 
ejecta does not deviate far from a spherical geometry, 
the relative masses and velocities of the two tails do no 
induce significant differences in the peak temperatures 



Electromagnetic Transients from NSNS/BHNS Mergers 



5 



of the two tails, and we have assumed a gray opacity. 
Therefore, we have shown that the uncertainty does not 
lie in the structure of the ejecta, but in the total ejected 
mass, its velocity structure, and the opacity of pure r- 
process material. The ejected mass and velocity are both 
depen dent on the nuclear EoS (Lee 2000; Oechslin et al. 



2007 ) and the treatment of gravity (Rosswog 2005), both 
of which we have not studied in detail here. Also, heating 
from r-process nucleosynthesis may effect the dynamics 



of the material and the total ejected mass (Metzger et al. 
|2010a|). The opacity of pure r-process material has not 
Deen studied at these low densities, so its relevance can- 
not be accurately assessed at this time. 

4. IMPLICATIONS FOR SGRBS, GW 
OBSERVATIONS, AND THE i?-PROCESS 

If compact object mergers are indeed the progenitors 
of SGRBs, these optical transients could possibly be seen 
along side their more spectacular prompt gamma-ray dis- 
play. In the second panel of figure |3j we compare our 
synthetic optical light curves to optical afterglow obser- 



vations and upper limits of SGRBs taken from Berger 
(|2010|). Clearly, the optical observations are all for on- 
axis SGRBs. We also show two synthetic on-axis SGRB 
optical afterglow l ight curves taken from )Van Eerten fc, 
MacFadyen (2011), and for comparison, the light curves 
of two standard supernovae. It seems that the optical 
observations can be reasonably explained by afterglow 
models, but, interestingly, they also do not rule out the 
possible contribution of an r-process powered supernova. 
It is also expected that in a sizeable fraction of events, 
these r-process powered supernovae may dominate the 
optical light at timescales of a day. 

The current and next generation of surveys (such as 
PTE, SASIR, Pan STARRS and LSST) wih be sensitive 
to very subtle changes in fiux, for sources with variability 
timescales of more than a few days, across most of the 
sky. In table [2j we show the detection rates expected 
for three transient surveys fflv en the low, r ecommended, 
and high merger rates given in | Abadie et_a l. (2010), cal- 
culated assuming a cosmology given by TComatsu et al.| 
(^2009) and a constant merger rate per volume. The rates 
are corrected for the effects of a finite cadence For the 
PTF and Pan STARRS, the detection rate is significantly 
lower than that expected for advanced LIGO whereas 
the inverse is true for SASIR and LSST. This suggests 
that, if these large numbers of events are observed, opti- 
cal observations will be more readily detectable than the 
associated gravity waves which could help focus gravita- 
tional wave searches and better constrain compact object 
merger rates jj Additionally, we find that LSST will po- 
tentially be able to detect these events out to z > 0.3, 
giving a handle on the evolution of the merger rate with 
cosmic time. 



It is clear from this and previ ous work ( Freiburghaus 



et al. 1999 Goriely et al. 2005) that large amounts of 



r-process material will be ejected in these events which 
may contribute significantly to the total r-process bud- 
get of the galaxy, although there may be some problems 

^ Even a cadence of 3 days reduces the detection rate by a factor 
of two compared to when the effects of a finite cadence are left out. 

^ We note that even if all compact mergers are associated with 
SGRBs, the r-process powered SNe will still dominate the off-axis 
optical afterglow. 



explaining the low metalli city halo r-process a bundance 
data within this scenario (Argast et al. 2004). Because 
the production of r-process elements is robust in the tails 
(i.e. the production of material with A > 120 is in- 
sensitive to the initial conditions of the material), the 
major uncertainty in this model lies in the amount of 
mass ejected in these events and h ow this material is 



spread through in dividual galaxies (Zemp et al. 2009^ 
Kelley et al.||2010 ). This is in contrast to the other pos- 



sible site of r-process nucleosynthesis, neutrino driven 
winds, where the uncertainty lies in finding conditions 
that make an r-process (c.f. Rober ts et al.||2010| ). For 
a given amount of mass ejection, the total amount of r- 
process material in the Milky Way puts a n upper limit 
on the merger rate (iMetzger et al. 2010b). Given these 
considerations, it seems likely that observations, or lack 
thereof of the transients produced in these mergers will 
give significant insight into the evolution of r-process ma- 
terial in the Universe. 
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Table 2 

Detection Rates for Various Blind Transient Searches. Properties of t he surveys are assumed to be the same as those given in |Strubbe 
Quataert| ( |2009| . 







PTF 






Pan STARRS 






LSST 






SASIR 


Model 


-^low 


Rre 


-^high 


-^low 


Rre 


-^high 


-^low 


Rre 


-^high 


-^low 


Rre Rhigh 


NS-NS 


0.1 


10 


100 


1.3 X 10" 


-3 0.13 


1.3 


2.3 X 10^ 


2.3 X 10^ 


2.3 X 10^ 


1.2 


120 1.2 X 10^ 


BH-NS 


7.3 X 10" 


-3 0.36 


12 


8.4 X 10" 


-5 4.2 X 10-3 


0.14 


1.6 


79 


2.6 X 10^ 


0.083 


4.1 140 



Note. — All rates are in detections per year. 
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